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ABSTRACT 



' Context. During an analysis of optical spectra of 80 young star clusters in several nearby spiral galaxies, [Om] and [Nn] emission lines were 
Xf \ . noted in some cases. Three of these emission line sources are associated with clusters older than 30 Myrs, and are identified as likely planetary 
■ nebula (PN) candidates. 

OO ' Aims. These objects may represent a rare opportunity to study PNe whose progenitor stars are known to be of intermediate masses, although 
detailed analysis is challenging because of the underlying strong continuum from the cluster stars. This paper presents and discusses basic 
properties of the PN candidates and their host clusters. 

Methods. Based on the observed emission line fluxes, the excitation parameters and luminosities of the nebulae are derived. This allows a 
(**] . crude placement of the central stars in two of the objects on the H-R diagram. Cluster ages and masses are estimated from broad-band colours 
O l' and by fitting model SSP spectra to the observed spectra. 

J. Results. The two PN candidates where central star luminosities and temperatures can be estimated are found to be consistent with post-AGB 
»| . model tracks for a central star mass of about 0.60 M Q . One of the host clusters has an age of 32-65 Myrs, corresponding to a main sequence 
' turn-off mass of M TO = 6.6-9.0 M Q . For the other cluster the age is 282^1-07 Myrs, corresponding to M T0 = 3.2-3.6 M . By estimating the 
£^ number of stars evolving off the main sequence per year, a total of 6 PNe are expected in the full sample of 80 clusters for a PN lifetime of 10 4 
• • ■ years. The factor of two disagreement with the actual observed number may be due, among other things, to uncertainties in PN lifetimes. It is 
interesting to note that all three PN candidates are associated with clusters which are more diffuse than average. 

Conclusions. While PNe have previously been found in some old globular clusters, the candidates identified here are among the first identified 
$H ' m y° un 8 star clusters. 

a ; 
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1. Introduction 

Star clusters and planetary nebulae (PNe) can both be observed 
at distances well beyond the Local Group, and are thus valuable 
extragalactic population tracers. However, much remains to be 
understood about the post-main sequence evolution of low- and 
intermediate stars leading up to the PN stage, partly because of 
the difficulty in determining the distances and progenitor stel- 
lar masses of planetary nebulae in the Milky Way. For PNe 
associated with stellar clusters the initial mass of the progen- 
itor star can be assumed to be that corresponding to the main 
sequence turn-off (unless the progenitor star is a close binary 

Send offprint requests to: S. S. Larsen, e-mail: 
larsen@astro .uu.nl 

* Based on observations collected at the European Southern 
Observatory, Chile under programme 072.B-0242(A) and with the 
NAS A/ES A Hubble Space Telescope, obtained at the Space Telescope 
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for Research in Astronomy, Inc., under NASA contract NAS 5-26555 



in which mass exchange has occurred), but the number of PNe 
known to be associated with stellar clusters remains small. Less 
than a handful of such objects are known to exist in Milky Way 
open clusters (Pedreros [TW7l[T^l 0'Dell[l963} and four PNe 
have been found in Galactic globular clusters (M15, M22, Pal 
6, NGC 6441; Jacoby et al. nWfl . 

The distance problem can be alleviated by studying PNe in 
external galaxies with well determined distances. Many such 
searches have been carried out or are currently under way 
(Kwok 2000, Feldmeier l2006l Magrini 2006 ), but the detected 
candidates are typically field objects and the initial mass thus 
remains unknown. PNe have been identified in a few GCs in the 
nearby giant elliptical NGC 5128 (Minniti & Reikuba 12001 
Rejkuba et al. 2003 1, and the significant number of on-going 
large spectroscopic surveys of extragalactic globular cluster 
systems are likely to reveal many more candidates. However, 
by studying PNe in ancient globular clusters one remains fun- 
damentally limited to low-mass stars. One way to sample a 
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larger range of stellar masses would be to collect a sample of 
PN candidates in younger stellar clusters. So far, the number of 
such objects are limited to the small number of PNe in Milky 
Way open clusters. 

We have recently obtained optical spectra for a number of 
young star clusters in several nearby galaxies. The original aim 
of this programme was to carry out a detailed test of simple 
stellar population (SSP) models, but during the analysis we 
noted emission lines in some of the spectra. As discussed be- 
low, the clusters in question span a range of ages from a few 
tens of Myrs to several hundreds of Myrs, corresponding to 
main sequence turn-off masses between 3 and 9 M . These ob- 
jects potentially represent a valuable addition to the very few 
cases where PNe have been associated with young star clusters, 
and could provide a unique opportunity to study PNe where 
the distances and initial masses are known with some accuracy. 
The full spectroscopic dataset and analysis will be presented in 
a forthcoming paper, and here we concentrate on the few spec- 
tra where PN candidates were identified. 

2. Data 

Spectra of a sample of young massive clusters (YMCs) in sev- 
eral nearby spiral galaxies were obtained with the EMMI in- 
strument on the ESO New Technology Telescope from March 
16 - March 18, 2004. The observations were carried out in 
multi-slit mode, using the RILD ("Red Imaging and Low- 
Dispersion Spectroscopy") mode and grism #3, which provided 
a spectral resolution of A/AA ~ 760 over the wavelength range 
4000A- 9000A. In each of the 7 slitmasks, typically 3 expo- 
sures of 3600 s each were obtained with the CCDs read out in 
2x2 binned mode. The full spectroscopic sample included 80 
clusters in four galaxies (NGC 2835, NGC 2997, NGC 3621 
and NGC 5236), selected primarily from the sample presented 
in Larsen ( 1999 1. In order to maximise the fraction of bona-fide 
star clusters (as opposed to "asterisms" or loose associations), 
as many objects as possible with imaging from the Hubble 
Space Telescope (HST) were included (Larsen 2004). The as- 
trometry for the slitmask design was based on pre-imaging data 
obtained with EMMI in service mode on Jan 28, 2004, exposed 
for 5 min in the R -band. 

The spectra were reduced with standard tools in the IRAF 1 
ONEDSPEC package and were wavelength- and flux calibrated 
using observations of HeAr calibration lamps and the flux stan- 
dards Hiltner 600 a nd LTT 2415 (Stone & Baldwin [T9831 
Hamuy et al. 1994). Small additional corrections (on order 
~ 1 A) to the wavelength scale were applied using the [Oi] 
5577. 35A night sky line. 

Accurate absolute flux calibration of slit spectra is always 
challenging due to slitlosses, which depend on the seeing and 
accurate centring of objects on the slit. Because EMMI does 
not have a functioning atmospheric dispersion corrector, an ad- 
ditional complication for multislit spectroscopy is differential 

1 IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc. under contract with the National 
Science Foundation 



refraction, which causes the slitlosses to be wavelength depen- 
dent as the slits cannot be aligned with the parallactic angle 
for the entire duration of an exposure. Other factors, such as 
the spatial distribution of the targets, also limit the freedom to 
choose the orientation of the slits. Although the observations 
were carried out at low airmass (generally < 1.2), some dif- 
ferences were noted in the flux calibration between individual 
observations of the same targets. 

We attempted to bring our spectra closer to an absolute 
flux scale by scaling them to match the fluxes derived from 
the ground-based photometry. Since the wavelength range cov- 
ered by the EMMI spectra includes four of the photometric 
passbands (BVRI) for which broad-band photometry is avail- 
able, the wavelength-dependent nature of the slit losses could 
also be accommodated. This was done by forming the ra- 
tio of the photometric fluxes (F p h ot (A)) in each band to the 
fluxes measured in the corresponding wavelength ranges of 
the spectra (F spectro (A)), and then producing a correction func- 
tion by fitting a second-order polynomial to the resulting ra- 
tio F p h ot /F spectro (A). Within a given exposure, differences in the 
overall scaling of individual spectra could be as large as ~ 25%. 
This may reflect uncertainties in the photometric as well as the 
spectroscopic measurements. The shapes of the individual cor- 
rection functions showed better agreement. Furthermore, for 
some objects (e.g. Hn regions dominated by strong emission 
lines) it was not possible to construct reliable correction func- 
tions. We therefore applied one average correction function to 
all spectra obtained in a given exposure, based on about half 
of the objects in a given slitmask. Typically, the correction 
amounted to 10%-20%, but in the most extreme cases it could 
be as much as a factor of 2. Each spectrum was then multiplied 
by the correction function, and the several individual spectra of 
each object were averaged (using a sigma-clipping algorithm 
to reject cosmic -ray hits) to produce the final spectra. 

Visual inspection of the reduced spectra generally revealed 
A-type spectra with strong Balmer lines, as expected for young 
star clusters with ages on the order of ~ 10 8 years. Some very 
young objects, known from ground-based narrow-band imag- 
ing to be embedded in Ha emission, were also included, and 
showed the expected strong emission lines from H, N, and 
O characteristic of Hn regions. Here we concentrate on three 
spectra which showed weaker [O in] and [N n] emission lines 
superimposed on A-type stellar spectra, for which both broad- 
band photometry and a comparison with model spectra suggest 
ages significantly greater than 10 7 years ( ^3.2> . Ground-based 
Ha images showed no line emission in the immediate vicinity 
of these objects, again suggesting that the emission-line objects 
contained in the star clusters are not H n regions. As discussed 
below, a more likely (although not entirely unique) explanation 
is that these clusters contain planetary nebulae (PNe). 

The three clusters discussed here are located in the two 
galaxies NGC 5236 and NGC 3621, both of which are ac- 
tively star-forming, nearby spiral galaxies. Throughout this pa- 
per, we assume distance moduli of m - M — 28.25 + 0.15 for 
NGC 5236 (T him et al. l2003l and 29.10 + 0.20 for NGC 3621 
(Rawson et al. 1997 1. For the foreground extinction we use the 
Schlegel et al. ffWSl values of A B = 0.284 mag (NGC 5236) 
and Ab = 0.346 mag (NGC 3621). For a more comprehensive 
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Fig. 1. EMMI spectra of the three YMCs with PN candidates. 
Each panel shows the raw spectrum (thick lines) and the residu- 
als when the best-fitting solar-metallicity model spectrum from 
Gonzalez-Delgado et al. (2005 1 is subtracted (blue lines). 

N5236-254 N5236-487 N3621-1106 




Fig. 2. HST/WFPC2 F606W (» V-band) images of the three 
YMCs. Each panel is 4" x 4" and the scale is 07 10 pixel -1 . 

listing of the properties of the host galaxies we refer to Larsen 
& Richtler 

3. Analysis 

3. 1 . Basic properties of the emission-line sources 

The EMMI spectra of the three objects are shown in Fig. ffl 
All three objects happened to have images from the Wide 
Field Planetary Camera 2 (WFPC2) on board the Hubble Space 
Telescope (HST), as shown in Fig. [5] Two of the objects are 
located in NGC 5236 (N5236-254, N5236-487 in the list of 
Larsen rr^ , while the third is in NGC 3621 (N3621-1 106). 
The morphological appearance of all three objects is generally 



Table 1. Emission line fluxes in units of 10~ 17 erg s -1 cm -2 . 
Note that both lines in the [O m] and [N n] doublets were fitted 
together. 





Fluxes 


[10 -17 erg s - 


cm- 2 ] 




N5236-254 


N5236-487 


N3621-1106 


Line 








H/3 4861 


6.0 


2.2 


4.6 


Ha 6562 


18.4 


12.3 


-0.15 


[Om] 4959,5007 


6.9, 20.6 


4.1, 12.1 


7.2, 21.7 


[Nn] 6548,6584 


5.6, 16.7 


0.7, 2.2 


4.7, 14.1 


Mean Error 


1.3 


0.8 


2.9 


"25007 


25.5 ± 0.06 


26.1 ±0.07 


25.4 ±0.13 



consistent with that of star clusters, all of them being signifi- 
cantly more extended than the WFPC2 point-spread function, 
but still fairly compact and with regular, symmetric profiles. 
Note, though, that N3621-1106 has a fainter companion at a 
projected separation of 077 (~ 22 pc). The companion is also 
resolved in the HST images, and thus also a likely cluster candi- 
date. Whether the two objects represent a chance alignment or 
a real physical pair is difficult to tell, however. Similar colours 
would suggest similar ages for the two objects and thus make 
a chance alignment less likely, but the ground-based imaging 
has insufficient resolution to address this question and the HST 
data are available in only one band. The half-light radii listed 
in Larsen J2004> are 10. 1 pc, 6.7 pc and 7.2 pc for N5236-254, 
N5236-487 and N3621-1106, respectively. While these sizes 
are consistent with those of star clusters, it is interesting to note 
that all three clusters have larger half-light radii than the mean 
value of 3-4 pc which is typical for both young and old clusters 
(Whitmore et al. rT999l Larsen l2004l Jordan et al. l2005l . 

The spectral signal-to-noise ratio per pixel in the disper- 
sion direction is about 106, 30 and 24 at 5000 A for N5236- 
254, N5236-487 and N3621-1106, respectively. The [Om] 
/14959,5007A and [Nn] ^6548,6584A lines are noticeable in 
all three spectra, but much weaker than typically seen in H n 
regions. Balmer lines are expected to be seen in emission in 
both planetary nebulae and H n regions, but in our cluster spec- 
tra the Balmer lines are dominated by the underlying strong 
absorption lines from the stellar component. 

In order to remove the stellar component from the spec- 
tra, model spectra from the library of Gonzalez-Delgado et al. 
( 2005 ) were subtracted from the observed spectra. The library 
contains model spectra for the integrated light of simple stellar 
populations with ages between 1 Myr and 17 Gyr, tabulated 
at 0.3A resolution. A Salpeter (1955)-like IMF is assumed. 
The best-fitting model spectra were selected by minimising the 
r.m.s. residuals between the observed spectra (shifted to zero 
radial velocity) and smoothed model spectra. Because of the 
uncertainties in the absolute flux calibration and the redden- 
ing corrections, the model spectra were multiplied by a poly- 
nomial to follow the same large-scale spectral energy distri- 
bution as the observed spectra. This was done by fitting the 
ratio of the observed to model spectra with a 5th-order poly- 
nomial and multiplying the model spectra with this polyno- 
mial, rather than scaling them by a single constant, before cal- 
culating the r.m.s. residuals. The Gonzalez-Delgado et al. li- 
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Fig. 3. Gaussian emission line fits to the spectra after subtrac- 
tion of best-fitting solar-metallicity SSP models. In each panel, 
the thick lines represent the observed spectra while the red 
lines are the Gaussian fits. The regions around HyS+[Om] and 
Ha+[Nn] are shown on the same flux scale. 



by fitting Gaussian line profiles to the emission line spec- 
tra. The [Ora] ^5007/4959A and [Nn] .16584/6548A line ra- 
tios were kept fixed at 3.0 as dictated by atomic physics (e.g. 
Osterbrock 1988 1, while the [O m]/HB and [N n]/Ha ratios were 
allowed to vary. The resulting best-fitting Gaussian models are 
overplotted on the spectra in Fig. [2] For the two emission-line 
sources associated with YMCs in NGC 5236, no further shifts 
in wavelength were required to match the line profiles and the 
sources are likely to be physically associated with the clusters. 
The best fit for N3621-1 106 was obtained if an additional shift 
of about +150 km/s was applied to the observed spectrum, so in 
this case the source of the emission lines may not be physically 
associated with the cluster. 

Table^hsts the emission line fluxes obtained from the fits 
in Fig. |3 Errors were estimated by fitting emission lines at 8 
different wavelengths where none were expected. The mean 
errors in the table represent the standard deviation of these 8 
test measurements. It should be emphasised that these error 
estimates only account for random measurement errors. The 
line fluxes may be further affected by systematic errors, for ex- 
ample due to inaccurate subtraction of the stellar continuum. 
To investigate this issue, we repeated the emission line mea- 
surements after subtracting Z = 0.008 SSP models based on 
the Padua isochrones, and Z = 0.008 and Z = 0.020 models 
based on the Geneva isochrones. The resulting emission line 
fluxes generally agreed with those listed in Tabled within the 
stated one-sigma errors. Only for the Hfi line in the spectrum 
of N5236-254 was a somewhat larger deviation encountered - 
for this line, a flux of 2.5 x 10~ 17 erg -1 cm" 2 was returned when 
using the Z = 0.008 Padua-based SSP model, i.e. a 2.7 'cr devia- 
tion. In conclusion, although some systematic errors cannot be 
entirely ruled out, the error estimates in Table appear to be 
fairly realistic overall. The flux calibration may account for a 
further ~ 10 - 20% uncertainty on the line fluxes, although rel- 
ative line ratios are probably more accurate. Finally, the msixn 
magnitudes, defined as OTgoo7 = -2.5 log(F5007)- 13.74 (Jacoby 
1989), are also listed. 



3.2. Host cluster properties 



brary includes models based on both the Geneva and Padua 
stellar isochrones and for a range of metallicities, but for this 
work we only used Solar metallicity models based on the Padua 
isochrones. The analysis carried out in the following does not 
depend strongly on the model choice. For NGC 5236, analysis 
of H n regions (Bresolin & Kennicutt 2002 1 yields an oxygen 
abundance very close to Solar [12+log(0/H)=8.9] at the galac- 
tocentric distance of the two clusters (6.'3 for N5236-254 and 
5.'1 for N5236-487). For NGC 3621, the work of Ryder ( fT995t 
indicates 12+log(O/H)=9.0 at the position of N3621-1106, V.O 
from the centre. The use of Solar metallicity models thus ap- 
pears justified. The only remaining free parameter is then the 
age, and for each cluster we selected the model spectrum of the 
age which gave the best fit to the data. 

Fig. [3] shows the regions around the Hj6+[Oin] and 
Hff+[Nn] lines after subtraction of the best-fitting, scaled 
solar-metallicity models. Emission line fluxes were measured 



Photometry for the clusters was taken from Larsen (1999 1. 
Table [2] lists the UBVRI broad-band colours for each clus- 
ter, corrected for foreground extinction only. From the broad- 
band colours, rough age estimates can be obtained by compar- 
ison with predictions by simple stellar population (SSP) mod- 
els, or by using the empirical 'S' -sequence calibration (based 
on LMC clusters) by Girardi et al. U9951 . Fig. |4] shows the 
U—B vs. B — V two-colour diagram, including the S-sequence 
and photometry for the three clusters. Based on this plot, 
N5236-487 appears to be the youngest of the three objects with 
log(age)~7.8, while N5236-254 is the oldest (log(age)- 8.6). 
Including all the colours available in Table |2] and comparing 
with solar metallicity SSP models from Bruzual & Chariot 
(2003 1, the photometric age and extinction estimates listed in 
Table |2 were obtained. Here, AAb denotes the extinction re- 
quired to best match the observed colours in addition to the 
foreground extinction. Errors on the age- and extinction esti- 
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Table 2. Cluster data. Photometry is corrected for foreground 
reddening only. AA B is the additional amount of reddening re- 
quired to best match the SSP model colours. The ranges of 
turn-off and total cluster masses reflect the differences between 
spectroscopic and photometric age estimates. The masses are 
for a Salpeter-like IMF with a lower mass limit of 0.1M Q 





N5236-254 


N5236-487 


N3621-1106 


V 


18.20 ±0.01 


19.96 ± 0.04 


19.67 ±0.14 


(U - B) 


0.16 + 0.02 


-0.33 ± 0.05 


-0.13 ±0.15 


(B - V) 


0.28 ± 0.01 


0.20 ± 0.05 


0.19 ±0.17 


(V-R) 


0.22 ± 0.01 


0.27 ± 0.05 


0.24 ±0.21 


(V-l)o 


0.46 ± 0.01 


0.60 ± 0.05 


0.57 ±0.23 


log(age/yr) phot 


8.56 + 0.03 


7.81+0.24 


8.21 ±0.45 


AA B [mag] 


0.00 ± 0.03 


0.45 + 0.29 


0.30 ±0.55 


log(age/yr) spec 


8.45 


7.50 


7.70 


M/10 3 M G 


376-495 


29-51 


94-187 


M T0 /M 


3.2-3.6 


6.6-9.0 


4.5-7.3 



-1.5 


6.55; 


I 




I 


-1 .0 




7.0 






m -0.5 




7.5 -A 










8.0 < 


L n5236-487 
n3621 -1 106 




0.0 




8 


.5 9.0 








n5236- 


254 • ' V ^S-^ 


9.5 


0.5 




i 




I 


-0.5 


0.0 


0.5 


1.0 



B-V 

Fig. 4. U—B vs. B— V two-colour diagram for the three YMCs. 
Also shown is the S-sequence from Girardi et al. (1995). The 
numbers along the S-curve indicate the log(age) corresponding 
to the major tick marks. 



mates were estimated by a Monte-Carlo procedure: First, ran- 
dom offsets (5 co i were added to the observed colours, and the 
age- and extinction were then re-derived. The offsets <J co i were 
drawn from a Gaussian distribution with standard deviation 
equal to the photometric errors on each colour. The errors on 
the ages and extinctions were then estimated as the standard de- 
viation of the individual values obtained from 100 such experi- 
ments. The ages derived from the comparison with SSP models 
agree quite well with those inferred from Fig.|4] 

The extinction estimates in Table [2] can be compared with 
those derived from the HafHB line ratios in Tabled Assuming 
the Case B recombination ratio F(Ha)/F(HB) = 2.85 (e.g. 
Osterbrock lT9"88l . we derive AA B = 0.02+0.95 mag forN5236- 
254 and AA B = 2.51 ± 1.53 mag for N5236-487. These values 
are consistent with the photometric AA B values, although the 



errors are larger. In the case of N3621-1 106, both Ha and WB 
are essentially non-detections, and it is not possible to derive a 
meaningful extinction estimate. 

Also listed in Table [2] are the spectroscopic age estimates, 
which are simply the ages of the best-fitting SSP model spec- 
tra. The spectroscopic ages tend to be lower than the photo- 
metric ones, although for any individual cluster the difference 
between the two age estimates is not much larger than the lcr 
error. However, the relative ranking of the clusters remains the 
same. A detailed comparison of photometric and spectroscopic 
ages, using the full sample, will be carried out in a separate pa- 
per. Because the continuum shape was not fitted, no constraints 
on the extinction are available from the spectral fits. 

Once the ages and reddenings are known, the cluster 
masses can be estimated by combining the SSP model predic- 
tions for mass-to-light ratios with the observed magnitudes and 
adopting the distances for the galaxies quoted above. Table |2] 
lists these photometric mass estimates, assuming a stellar mass 
function with the Salpeter slope {dN/dM oc M~ 235 ) extend- 
ing down to 0.1 M . For a more realistic IMF (e.g. Rroupa 
2002 1 the masses would be lower by a factor of ~ 1 .6. The mass 
ranges given in the table correspond to the difference between 
the photometric and spectroscopic ages, with the lower masses 
corresponding to the younger ages. 

The main sequence turn-off masses, Mto, corresponding to 
the cluster ages were obtained from the Salasnich et al. (2000 1 
stellar model tracks. Again, Table |2] lists a range correspond- 
ing to the photometric and spectroscopic ages. If the emission 
line sources are indeed PNe associated with their putative host 
clusters, then their progenitor stars should have started out with 
masses close to the cluster turn-off masses. 

3.3. The nature of the sources: planetary nebulae? 

Several classes of objects show emission line spectra, but not 
all are consistent with the characteristics of the objects studied 
here. First, H n regions are common in star-forming galaxies. 
Our full sample contains several such objects, which generally 
show very different spectra than those in Fig.^ The emission 
line equivalent widths of H n regions are generally much larger 
than in the spectra in Fig. [2 and the excitation is typically lower 
(e.g. Acker et al. 1987 >. Furthermore, young star-forming re- 
gions tend to have a much less regular structure than the smooth 
profiles indicated by the images in Fig.|2] Finally, only clusters 
younger than ~ 10 Myrs contain stars which are hot enough to 
produce the ionising radiation required to excite an H n region, 
while we derive significantly higher ages for the three objects 
studied here C ^3.2> . We can therefore exclude with fairly high 
confidence that the objects are H n regions. 

Wolf-Rayet stars also show strong emission lines, but again 
the spectra are quantitatively (and qualitatively) very different 
from those observed here. In particular, the lines in W-R star 
spectra are very broad (several 1000 km/s) and would be eas- 
ily resolvable at the spectral resolution of our EMMI data. As 
for H n regions, the ages of the clusters studied here make it 
highly unlikely that they contain W-R stars, whose progenitors 
are believed to be massive (> 2OM ; Humphreys et al. 1985 1 
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stars with lifetimes of only a few Myrs. We note, however, that 
a few objects in our full sample do show the characteristic W-R 
bumps at 4650A and 5812A (e.g. Torres & Massev [T3g7l . but 
these features are entirely absent from the spectra in Fig.^ 

A third possibility, which is more difficult to dismiss, is 
symbiotic stars (SSs). These are binary systems in which a red 
giant is transferring material to a hot compact companion, such 
as a white dwarf (Kenyon 1986 1. Mass loss from the red gi- 
ant results in a nebula which is ionised by the hot companion, 
producing a spectrum very similar to that of a planetary neb- 
ula. The main difference is that the spectrum of a SS contains 
an extra component in addition to the emission-line spectrum 
and the spectrum of the hot star, namely that of the red giant. 
However, for our sources it is impossible to detect this small 
addition to the continuum on top of the cluster spectrum. On 
the basis of the data at hand, it is difficult to distinguish be- 
tween symbiotic stars and planetary nebulae as the most likely 
sources of the emission lines seen in our spectra. The progeni- 
tor stars are also expected to be similar, in both cases involving 
red giants or AGB stars. 

Most surveys of Galactic and extra-galactic emission-like 
objects identify a larger number of PNe than SSs, and we might 
thus expect on statistical grounds that our candidates are more 
likely to be PNe. About 1500 PNe are known in the Milky 
Way (Kohoutek 2001 1, while less than 200 symbiotic stars are 
catalogued (Belcyhski et al. 12000ft . Unfortunately, the com- 
pleteness of these catalogues is very difficult to assess for ei- 
ther type of object, and the classification of at least 50 objects 
in the Milky Way alone is still ambiguous (Kohoutek 2001 1. 
The total number of PNe in the Milky Way is estimated to be 
about 25000 (Buzzoni et al. 2006), but this is little better than 
an order-of-magnitude estimate. The uncertainty on the num- 
ber of SSs may be even greater, with estimates ranging from 
3000 rKenvon [T986T l to as many as 300.000 (Munari & Renzini 
ri992l 

Surveys of PNe and SSs in the Magellanic Clouds have the 
advantage that the distances are known with reasonable accu- 
racy. About a dozen SSs are known in the LMC and SMC. 
These have H/3 luminosities in the range 7.7 x 10 33 - 5 x 10 34 
erg s~' (Morgan 1992 1, i.e. all are fainter than the H/3 fluxes 
derived for the objects studied here ( 33.41 1. For typical effec- 
tive temperatures of the hot component of ~ 10 5 K (Miirset et 
al. fT996l . the [Om] ,15007 flux is about an order of magnitude 
greater than the H/3 flux. As a crude estimate, there may then 
be about 10 SSs in the Clouds with L([Om]) > 10 35 erg s _1 
(although, again, an unknown completeness correction may ap- 
ply). For comparison, more than 100 PNe with L(fO in]) > 10 35 
erg s _1 are catalogued in the Clouds (Jacoby et al. 1199 1. Thus, 
at least at the high-luminosity end, PNe appear to outnumber 
SSs by about an order of magnitude in the Magellanic Clouds. 
Furthermore, Miirset et al. & 1996ft conclude that the LMC SSs 
most likely belong to a population older than ~4 Gyrs, which is 
much older than the clusters studied here ( 33.2ft . It is also worth 
noting that no SSs are known to be associated with stellar clus- 
ters in the Milky Way while about 7 PNe have been found in 
Galactic globular or open clusters, as noted in the Introduction. 

We thus proceed under the assumption that the emission- 
line objects found here are most likely planetary nebulae, al- 
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log T eff [K] 

Fig. 5. H-R diagram for the PN candidates in NGC 5236, with 
log T e ff and log L cs /L estimated from the [Om] and H/3 lines 
using the calibrations given by Dopita et al. (1992). Solar- 
metallicity post-AGB hydrogen-burning evolutionary tracks 
from Vassiliadis & Wood (1994) are shown for comparison, 
with the main sequence and core mass (in parentheses) at 10 4 
K indicated. Both objects are consistent with a core mass close 
to 0.60 M . 



though symbiotic stars cannot be entirely ruled out as an al- 
ternative. Repeat observations might help clarify this issue, as 
SSs often show significant variability on time scales of days, 
weeks or months. A detailed analysis of the nebular properties 
using high-dispersion spectroscopy might also help distinguish 
between SSs and PNe, since SSs typically have ~ 100 times 
higher densities than PNe (Kenyon 1986). As a further note, it 
seems to us that the identification issue frequently receives less 
attention than it merits in extragalactic PN surveys. 

3.4. Physical properties of the PN candidates 

Because only the strongest emission lines can be measured in 
the spectra, it is not possible to derive basic nebular proper- 
ties such as temperature and density directly. While this lim- 
its the amount of physical information that can be extracted 
from the spectra, some constraints can be put on the total lu- 
minosities and effective temperatures of the ionising sources. 
Dopita et al. (1992) have published an extensive set of model 
calculations, allowing estimates of these properties from ob- 
servations of the Balmer lines and strong forbidden lines. They 
define the excitation parameter as E = 0.045 {F (5001) / F( Hp)} 
for 0.0 < E < 5.0. For higher excitation parameters the Hen 
/14686A line is required. This line is not detected in our spec- 
tra, but as will be seen below the excitation parameters derived 
here are sufficiently low that only the O and H lines are needed. 

Measurements of the Balmer emission lines rely critically 
on correct subtraction of the underlying cluster spectra. From 
Fig-E]it is difficult to assess the reliability of this subtraction for 
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Table 3. Derived nebular properties. L cs is the luminosity of the central star. Note that H/3 fluxes are derived from the Ha fluxes 
assuming standard case B recombination. 





N5236-254 


N5236-487 


N3621-1106 


L(Hor) [10 35 ergs 1 ] 


5.14±0.85 


4.44+0.73 




L(H0) [10 35 ergs" 1 ]* 


1.80+0.30 


1.56+0.25 




L ([0 III]) 


6.11+0.99 


5.10+0.84 


18.7+4.49 


M5007 


-3.01+0.18 


-2.82 + 0.18 


-4.23 ± 0.26 


Excitation parameter (E) 


1.52+0.14 


1.47+0.14 




lo g r eff [K] 


4.66+0.02 


4.65+0.02 




log L cs [L Q ] 


3.79+0.08 


3.73+0.08 





H/3. The Ha line, on the other hand, is stronger and less diluted 
by the stellar absorption, and therefore easier to measure reli- 
ably, at least for the two objects in NGC 5236. For the object in 
NGC 3621, neither H/3 nor Ha can be reliably measured. Thus, 
in the following discussion we concentrate on the two objects 
in NGC 5236, for which the association with the star clusters 
is anyway more secure. 

Table |3] lists the basic parameters derived for the nebulae. 
The H/3 fluxes were estimated from the measured Ha values as- 
suming the Case B recombination ratio, but the results do not 
change much if we use the measured H/3 fluxes directly. For 
N5236-254, the de-reddened H/3 flux [(1.80 + 0.47) x 10 35 erg 
s -1 ] agrees almost exactly with that derived from the Ha flux. 
For N5236-487 the de-reddened H/3 flux is (0.95 + 0.37) x 10 35 
erg s _1 , about 40% lower than the value derived from Ha, 
but with a larger formal error. All line fluxes have been de- 
reddened using the Cardelli et al. ( 1989 1 extinction curve and 
the extinction estimates from Table 13 The quoted error esti- 
mates include a 0.15 (0.20) mag uncertainty on the distance 
modulus of NGC 5236 (NGC 3621). For reference, the [Oin] 
luminosities have been converted to absolute M5007 magni- 
tudes. These are well within the range expected for planetary 
nebulae; the tip of the planetary nebula luminosity function is 
at M5007 ~ -4.5 (e.g. Kwok l2o6ol 

Based on the excitation parameters and the Hp luminosities, 
the effective temperatures and total luminosities of the central 
stars (L cs ) were derived using the relations given in Dopita et 
al. ([[9921. The errors on T e ff and L cs were estimated simply 
by propagating the random measurement errors, and thus do 
not include any systematic errors in the calibration. Such errros 
might result, for example, from metallicity differences between 
the LMC and SMC calibrator PNe used by Dopita et al. and the 
objects studied here. For the relevant temperature range, Dopita 
et al. find that the L(Hy8)/L cs ratio changes by only ~ 1 percent 
for O abundances between 0.1 and 2 times Solar (their Fig. 1). 
Over the same metallicity range, log T e g varies by ~ 0.05 dex 
for a fixed E S= 2. The error estimates in Table[5]do not include 
uncertainties on the reddening. For E and T e s the effect on the 
error budget is negligible (Alogr e ff = 0.01 mag for N5236- 
487) while the additional uncertainty on L cs is A log L cs = 0.06 
dex for N5236-487. In both cases, the associated errors for 
N5236-254 are even smaller. 

The H-R diagram for the two nebulae in NGC 5236 is 
shown in Fig. |5] together with H-burning post-AGB evolution- 
ary tracks from Vassiliadis & Wood ( 1994 1. Curiously, the two 
central stars have very similar luminosities and effective tem- 



peratures, in both cases corresponding to a mass close to 0.60 
M . This is somewhat lower than predicted by the Vassiliadis & 
Wood tracks for initial masses of 3.2-3.6 M Q (for N5236-254) 
and > 6.6M (N5236-487). Empirically, the initial-final mass 
relation for AGB and post-AGB evolution shows a fairly large 
scatter. It may be nearly flat for Mi n ; t i a i < 3 - 4 M , although 
most studies suggest that a star with an initial mass > 6.6M is 
more likely to produce a ~ 0.8 -1.0 M remnant (Kwok 2000 
Ferrario et al. 2005 Weidemann 20001. It would clearly be de- 
sirable to put stronger constraints on the physical properties of 
these PNe. 

We have used the solar-metallicity (Z=0.016) tracks in 
Fig. [5] but the Z=0.004 tracks do not differ much and the in- 
ferred core masses would be lower by only ~ 0.02 M . The 
formal errors from Table [3] are smaller than the plot symbols, 
but we emphasise that a more robust analysis, using better data, 
would be desirable. However, since the model tracks are nearly 
horizontal in the H-R diagram, uncertainties on T e g are of mi- 
nor importance for the comparison with the post-AGB tracks. 
Likewise, uncertainties on the flux calibration and reddening 
correction should lead to shifts of no more than 10-20%, or 
-0.1 dex in logL cs . The difference between the measured and 
derived H/3 flux for N5236-487 corresponds to a shift of 0.2 
dex (downward) in log L cs . 

From Fig. [3] one notes significant differences in the rela- 
tive strengths of the [Nn] and [Om] lines. In N5236-254 and 
N3621-1106 the [Nn] lines are only slightly weaker than the 
[Om] lines, while in N5236-487 they are about factor of 6 
weaker. This may appear somewhat puzzling, considering the 
very similar parameters derived for the central stars in the two 
PNe in NGC 5236. However, the PNe studied here fall well 
within the range of [O m]/[N 11] line ratios seen in Galactic PNe 
(Cuisinier 1996 1, which indeed show a very large range in the 
relative abundances of N and O. This is likely related to CNO 
processing in the progenitor stars (Peimbert & Serrano 1980 1 
and the ratio of N to O is expected on theoretical grounds to 
increase steeply with progenitor mass (Renzini & Voli 1981 1. 
With the current data we cannot carry out an actual abundance 
analysis, but it may be worth noting that the spectrum having 
the weakest N lines is that associated with the youngest cluster 
(and thus the highest progenitor star mass). The Ha line is also 
stronger relative to the [N 11] doublet in N5236-487 than in the 
other objects. At least qualitatively, this seems consistent with 
the overall decrease in the [N/O] abundance ratio with decreas- 
ing [Nn]/Ha line ratio noted by Perinotto & Corradi ( 1998 1, 
again hinting that N5236-487 may have a decreased [N/O] ra- 
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Fig. 6. Absolute My magnitude versus log(age) for all clusters 
in the spectroscopic sample. The clusters containing PN can- 
didates are indicated with filled circles. Ages in this figure are 
estimated from broad-band colours. 




Fig. 7. Stellar lifetime as a function of initial mass, according to 
Salasnich et al. (2000 1 model tracks. The curve is a power-law 
fit to the points. 



tio, rather than an enhanced one as expected for a massive pro- 
genitor. 

3.5. Observed vs. predicted number of PNe 

Fig. [6] shows the absolute magnitude My vs. the (photomet- 
ric) age estimate for all 80 clusters included in the full spectro- 
scopic sample. Ages were estimated by matching the observed 
UBVRI colours against SSP model predictions from Bruzual 
& Chariot d^UOll as described in Of the 80 clusters, 28 
have ages less than log(t/yr) = 7.5 and are thus unlikely to host 
planetary nebulae. For the remaining 52 clusters, it is of interest 
to compare the expected number of PNe with the 2-3 objects 
actually detected (in the clusters marked by filled circles). 




8.5 9.0 
log oge [yr] 

Fig. 8. Predicted number of planetary nebulae per 10 6 M ver- 
sus age for a single-aged stellar population, assuming a PN 
lifetime of 10 4 years, a Salpeter-like IMF, and the stellar mass- 
lifetime relation from Fig. [7] 
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Fig. 9. Predicted number of planetary nebulae per 10 6 L Q ,y ver- 
sus age for a single-aged stellar population, using the same as- 
sumptions as in Fig.|H]and mass-to-light ratios from Bruzual & 
Chariot SSP models. 



Recently, Buzzoni et al. ( 2006 1 have given an extensive dis- 
cussion of planetary nebulae as population tracers. However, 
their main focus is on populations older than about 1 Gyr and 
their results are not directly applicable to our case. In the fol- 
lowing it is assumed that a stellar cluster consists of stars born 
in a single burst of negligible duration compared to the age of 
the cluster. As the age spread within a cluster is typically only 
a few Myrs at most, this is a reasonable assumption for clusters 
which are old enough to host PNe. If the (initial) mass of a TP- 
AGB star which is just about to become a PN is labelled Mp, 
then the stellar death rate is the number of stars with masses in 
a small range (Mp..Mp + dM), divided by the difference in the 
stellar lifetime (t*) across the mass interval dM. The predicted 
number of PNe in the cluster is then 



^VpN.pred. - /PNTPN 



N*(Mp + dM)-N*(Mp) 
T*(M P ) - T*(M P + dM) ' 



(1) 
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where / PN is the fraction of stars which will actually end their 
lives as PNe and r PN is the observable PN lifetime. N*(M) is 
the number of stars in the cluster with masses less than M and 
t*(M) is the lifetime of a star with mass M. Eq. Q can be 
rewritten as 



M>N,pred. - /pN TpN 



M^(M P ) 



(dT*/dM)(M P ) 



(2) 



where £ is the stellar IMF, £(M) = dN* /dM, (normalised to 
a total mass of 1) and M c \ is the total mass of the cluster (see 
also Renzini & Buzzoni ll 9861 and Buzzoni et al. 2006 1. For a 
Salpeter ( 1955 1 IMF populated with masses between M m [ n and 



a + 2 



f a+2 



ja+2 



M" 



(3) 



with a = -2.35. The denominator in Eq. can be estimated 
from stellar evolutionary tracks. Fig. shows the stellar life- 
time until the TP-AGB phase as a function of initial mass ac- 
cording to the evolutionary tracks in Salasnich et al. (2000 1. 
Also shown in the figure is a power-law fit to t* vs. M, which 
gives t* = 1.00 x 1O 1O (M/M )~ 2 - 75 yr. Inserting the derivative 
of this fit in Eq. 0, the curve in Fig.[8]is obtained. This shows 
the predicted number of PNe per 10 6 M as a function of age, 
assuming that all stars go through the PN phase (/pn = 1) and 
that the PNe are observable for a period of tpn = 10 4 years 
(e.g. Marigo et al. l2001> . independent of initial mass. The life- 
time tpn is probably the largest contributor to the uncertainty 
in this calculation. We also note that the power-law fit tends to 
underestimate the slope at the low-mass end, so that the num- 
ber of PNe is over-predicted for ages greater than several Gyrs. 
However, for the age distribution of our clusters this is not a 
problem. 

Using mass-to-light ratios from the Bruzual & Chariot SSP 
models, the number of PNe per unit luminosity can be calcu- 
lated. This is shown in Fig. [9] for the V-band. While the ex- 
pected frequency of PNe still depends on age, the dependence 
on assumptions about the low-mass end of the IMF disappears 
when normalising to luminosity instead of mass. This is be- 
cause the PNe progenitors are among the most luminous stars 
in a SSP while the low-mass stars contribute very little to the lu- 
minosity of a SSP. This makes the curve in Fig.|9]more directly 
comparable to observations. Of course, the same IMF has to be 
used when computing the number of PNe per mass (Eq.|2j and 
for the subsequent evaluation of the mass-to-light ratios. Here 
we are assuming a Salpeter-like IMF extending down to 0.1 
M . The PN formation rate derived here agrees within about a 
factor of two with the estimates by Renzini & Buzzoni ( 1986 1. 

We can now estimate the number of PNe expected in each 
cluster in Fig.|6] From Fig. [8] and Fig.|9]it is clear that, on av- 
erage, less than 1 PN is expected per cluster. By summing up 
the expected numbers of PNe in all clusters in the sample we 
estimate a total of 6 PNe in the age range 7.5 < log(age) < 9.0. 
This is to be compared with the 2-3 objects actually observed. 
Note that this assumes an observable PN lifetime of 10 4 years 
- if, for example, the lifetimes are shorter by a factor of two, 
then the predicted number of PNe would be lower by the same 



factor. Since the PN lifetimes are a steep function of the cen- 
tral star mass, with more massive progenitor stars expected to 
produce shorter lifetimes (Schonberner & Blocker 1996 1, the 
agreement appears satisfactory and these rough estimates cer- 
tainly suggest that the detection of a small number of PNe in 
a spectroscopic survey like the one carried out here is indeed 
expected. Alternatively, we can compare with the empirical es- 



timate of the PN specific frequency from Magrini et al. ( 2003 1, 
who estimate an average of one PN per 10 6 92 V-band L for 
Local Group galaxies (ignoring any age dependence). The in- 
tegrated V-band luminosity of all clusters in our sample older 
than 30 Myrs is about 3.3 x 10 7 L , corresponding to 4 PNe. A 
search of PNe in a larger sample of young and intermediate-age 
star clusters might prove rewarding. 



4. Summary and conclusions 

In a sample of optical spectra of 80 young extragalactic star 
clusters, 3 emission line objects were detected which we iden- 
tify as likely PN candidates. Two of these have radial velocities 
consistent with those of their host clusters, while the third dis- 
plays a shift of about 150 km/s and may be a field object falling 
on the slit by chance. The analysis of the candidate PN spec- 
tra is complicated by the relatively low spectral resolution and 
the strong underlying stellar spectra. This is a problem espe- 
cially for the Balmer lines, which are very strong in absorption 
in the cluster spectra and thus mask the comparatively weak 
line emission from the PNe. However, it was found that the 
line emission can be reasonably well separated by subtracting 
out synthetic SSP model spectra from the library of Gonzalez- 



Delgado(2005l. 

All three objects have [Oin] emission line fluxes con- 
sistent with those expected from PNe. For the two objects 
where the association with clusters is most robust (both in 
NGC 5236), the excitation parameters and total luminosities 
were estimated and the central stars were placed on the H-R di- 
agram. Comparison with post-AGB evolutionary tracks from 
Vassiliadis & Wood ( 1994 1 indicates central star masses of 
about 0.60 M Q . 

Based on stellar evolutionary timescales, we estimate that 
a total of about 6 PNe are expected in the 80 star clusters for a 
PN lifetime of 10 4 years. Instead, 2-3 are detected here, which 
we consider reasonable agreement. Scaling the empirical esti- 
mate of the PN specific frequency from Magrini et al. ([2003 1 
predicts a total of 4 PNe in our sample, although the age dis- 
tribution of our PN host clusters is likely different from that of 
the Local Group stellar populations on which the Magrini et al. 
estimate is based. A systematic search for PNe in a larger sam- 
ple of extragalactic star clusters might put useful constraints 
on the PN lifetimes and the range of stellar masses that pro- 
duce PNe. High-dispersion spectroscopy would allow a better 
subtraction of the underlying cluster spectra and thus more ac- 
curate measurements of the emission lines, and would provide 
constraints on the kinematic properties (expansion velocities) 
of the nebulae. Further observations would also provide con- 
straints on time variability, and thus help distinguish between 
PNe and other possibilities (e.g. symbiotic stars). 
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